Liquid stability in a model for ortho-terphenyl 
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We report an extensive study of the phase diagram of a simple model for ortho-terphenyl, focusing 
on the limits of stability of the liquid state. Reported data extend previous studies of the same 
model to both lower and higher densities and to higher temperatures. We estimate the location of 
the homogeneous liquid-gas nucleation line and of the spinodal locus. Within the potential energy 
landscape formalism, we calculate the distributions of depth, number, and shape of the potential 
energy minima and show that the statistical properties of the landscape are consistent with a 
Gaussian distribution of minima over a wide range of volumes. We report the volume dependence 
of the parameters entering in the Gaussian distribution (amplitude, average energy, variance). We 
finally evaluate the locus where the configurational entropy vanishes, the so-called Kauzmann line, 
and discuss the relative location of the spinodal and Kauzmann loci. 



I. INTRODUCTION 

Recent years have seen a strong development of numer- 
ical and theoretical studies of simple liquid models, at- 
tempting to develop a thermodynamic description based 
on formalisms which could be extended to deal also with 
out-of-equilibrium (glassy) states QSSSUHSEIi. 
Hard spheres, soft spheres, Lennard Jones mixtures and 
simple molecular liquids E3, E3, O 113, EM EH E3l 
have been extensively studied. Within the potential en- 
ergy landscape (PEL) E3 thermodynamic approach, de- 
tailed comparisons between numerical data and theoret- 
ical predictions have been performed. Estimates of the 
number fl of local minima -basins- as a function of the 
basin depth and of their shape have been recently evalu- 
ated for a few models [IJ El El El El SllM M E! 
and from the analysis of experimental data [23, |2J, |25| . 
The PEL approach, which is particularly well suited for 
describing supercooled liquids, provides a controlled way 
to extrapolate the thermodynamic properties of the liq- 
uid state below the lowest temperature at which equili- 
brated data can be collected. For example, estimates of 
the locus where the configurational entropy vanishes, the 
so-called Kauzmann locus, can be given 26]; the Kauz- 
mann locus provides a theoretical limit to the metastable 
liquid state at low temperatures. Another limit to the liq- 
uid state is met on superheating and stretching the liq- 
uid, when the nucleation of the gas phase takes place. In 
this case a convenient way to define the limit of stability 
of the liquid state against gas nucleation is provided by 
the spinodal line, i.e.. the locus of point where the com- 
pressibility diverges [23. The Kauzmann line and the 
spinodal line define the region of phase space where the 
liquid can exist in stable or metastable thermodynamic 
equilibrium. 

Recent theoretical and numerical work has focused on 
the thermodynamic relation between these two curves m, 
I2H l29t I30L l3l| . A recent thermodynamic analysis |30ll31| 
suggests that the spinodal and the Kauzmann loci meet 
-in the (P, T) plane- with the same slope at a point cor- 
responding to the maximum tension that the supercooled 



liquid can sustain. Simple models which can be solved 
analytically, as hard E3 or soft El H3 spheres com- 
plemented by a mean field attractive potential, support 
such prediction. A numerical study of a Lennard Jones 
mixture is also consistent p9| . 

In this paper we consider the Lewis and Wahnstrom 
rigid model for the fragile glass former ortho-terphenyl 
(OTP) [13, whose dynamic [H E! E3 and thermody- 
namic features [13, H3 have been studied in detail. Our 
aim is to calculate the spinodal and the Kauzmann lines 
to estimate the region of stability of the liquid, and study 
the relation between these two loci. 

We improve the data base of phase state points 
of previous studies |17| . extending to both lower and 
higher densities and to higher temperatures. Perform- 
ing analysis of the pressure- volume relation along several 
isotherms, we estimate the homogeneous nucleation line 
and the spinodal curve; we also report upgraded esti- 
mates of the statistical properties of the landscape sam- 
pled by the liquid. We confirm that at all densities a 
Gaussian landscape properly models the thermodynam- 
ics of the system in the supercooled state. Such agree- 
ment gives us confidence in the evaluation of the locus 
along which the configurational entropy vanishes. We fi- 
nally discuss the limits and possibilities of an analysis 
based on the inherent structures thermodynamic formal- 
ism. The description of our results is preceded by a short 
review of the potential energy landscape approach to the 
thermodynamics of supercooled liquids. 

II. BACKGROUND: THE FREE ENERGY IN 
THE INHERENT STRUCTURES 
THERMODYNAMIC FORMALISM 

An expression for the liquid free energy in the range 
of temperatures where the liquid is supercooled can be 
given within the PEL formalism. In supercooled states, 
i.e., when the correlation functions show the two steps 
behavior typical of the cage effect |36|, |3j| , the system's 
properties are controlled by the statistical properties of 
the PEL [33 ■ In the PEL formalism, the potential energy 
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hyper-surface -fixed at constant volume- is partitioned 
into basins; each basin is defined as the set of points 
such that a steepest descent path originating from them 
ends in the same local minimum. The configuration cor- 
responding to the minimum is called inherent structure 
(IS), of energy eis and pressure Pis- The partition func- 
tion can be expressed as the sum over all the basins, 
weighted by the appropriate Boltzmann factor, i.e., as 
a sum over the single basin's partition functions. As a 
result, the Hclmholtz liquid free energy F(T, V), at tem- 
perature T and volume V, can be written as |18j : 



terized using the expression: 

{6N-3 ~| 
^ HphuJiiEjs)) W = a(V)+b(V)E IS +c(V)El 

_(3) 

The evaluation of f an h(T, V) is described in Ref. 43] . 

Finally, S con f can be calculated as the difference of 
the (total) entropic part of Eq. and the vibrational 
contribution to the entropy: 

(T, V) = S(T, V) - 5 harm (T, V) - S anh (T, V). (4) 



F(T, V) = E IS {T, V) - TS conf (T, V) + f vib {T, V). (1) 



Here, Eis is the average energy of the IS explored at 
the given (T, V), f V ib is the vibrational free energy, i.e., 
the average free energy of the system when constrained 
in a basin of depth eis, and S con f(T, V) is the config- 
urational entropy, that counts the number of explored 
basins. S con f(T,V) is a quantity of crucial interest, 
both for comparing numerical results with the recent 
theoretical calculations 0, |3!| , and to examine some of 
the proposed relations between dynamics and thermo- 
dynamics 0, EH 03- Therefore, in order to evaluate 
the free energy one needs to estimate the three terms of 
Eq. • Eis (T, V) is calculated by means of a steepest 
descent potential energy local minimization of equilib- 
rium configurations (see Ref. for details). In frag- 
ile liquids the T dependence of Eis {T, V) follows a 1/T 

law q mm m. 

The basin free energy f v u, (T, V) takes into account 
both the basin's shape -curvature- and the system ki- 
netic energy. From the formal point of view, this term 
is the integral of the Boltzmann factor constrained in 
a basin, averaged over all the basins with same depth 
eis(T, V). Numerically, f v n, is evaluated as sum of two 
contributions: i) a harmonic contribution, which depends 
on the curvature of the accessed basins corresponding to 
the minimum at the given eis', w) an anharmonic con- 
tribution, usually approximated as a function of T only. 
In the case of a rigid molecule model, the harmonic con- 
tribution is the free energy associated with a system of 
(6N — 3) independent oscillators of frequency Wk, where 
the ojk are the square root of the eigenvalues of the sys- 
tem Hessian matrix -the matrix of the second derivatives 
of the potential energy- evaluated at the basin minimum 
(see R ef. [431 for details). This contribution can be writ- 
ten as |45j: 



{6JV-3 ~| 
HP^ieis)) W, 

(2) 

where the symbol ()' denotes the average over all the 
basins with the same energy eis- The Eis dependence 
of fvibi in the harmonic approximation, can be parame- 



III. BACKGROUND: THE GAUSSIAN 
LANDSCAPE 

In order to evaluate analitically the free energy in the 
PEL formalism, it is necessary to provide a model for 
the probability distribution of eis , i.e., the number 
Q{fiis)deis of basins whose depth lies between eis an d 
£is + deis- Among several possibilities, the Gaussian dis- 
tribution [2l|, |43, H3 seems to provide a satisfactory de- 
scription of the numerical simulations of Refs. |li . ll7ll47|] . 
The "Gaussian landscape" is defined by: 
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-deis- 



(5) 



Here, the amplitude e aN accounts for the total number 
of basins, E Q plays the role of energy scale and a 2 , mea- 
sures the width of the distribution. One can grasp the 
origin of such distribution invoking the central limit the- 
orem. Indeed, in the absence of a diverging correlation 
length, in the thermodynamic limit, each IS can be de- 
composed in a sum of independent subsystems |4Sj , each 
of them characterized by its own value of eis ■ The IS en- 
ergy of the entire system, in this case, will be distributed 
according to Eq. (JSJ. 

The assumptions of a Gaussian Landscape (Eq. ©) 
and of a quadratic dependence of the basin free energy 
on Eis (Eq. ©) fully specify the statistical properties of 
the model. Thus, it is possible to evaluate the T depen- 
dence of Eis and S con f. The corresponding expressions 
are [13 : 

p rrT/ s (E (V)~(b(V) + f3)gZ(V)) B 

Els[T ' v) = lTW) = T' 

2 (6) 

where, for convenience, we have defined A = "f"" 6 ^ and 

7 ' l+2ccH 



B 



and 



Sconf (T,V)/k B = a{V)N 2a 2 (V) ' ^ ' 

Note that Ei S is linear in 1/T. The predicted 1/T de- 
pendence of Eis and the parabolic dependence of S con f 
has been confirmed in several models for fragile liq- 
uids E1II3- 



FIG. 1: Results of the landscape analysis for the densities 
198 and 200 cm 3 /mol (squares and circles symbols respec- 
tively): a) average inherent structure energy (symbols) to- 
gether with the 1/T dependence of Eq. @ (solid line), b) 
anharmonic contribution to the entropy, c) harmonic free en- 
ergy ln(exp I^D^ -3 ln(/3ftaJi(_E/s))|), d) configurational en- 
tropy S CO nf together with the fits according to Eq. Q. 



From the relations above and fits of the numerical data, 
one obtains: i) the vibrational coefficients a, b and c from 
Eq. @; ii) the distribution parameters, E Q and <J 2 (V) 
from Eq. ©; Hi) the amplitude e aN from Eq. Q. A 
study of the volume dependence of the parameters a(V), 
E Q (V), and <J 2 (V) , associated with the F-dependence of 
the shape indicators (Eq. ©), provides a full character- 
ization of the volume dependence of the landscape prop- 
erties of a model, and offers the possibility of developing 
an equation of state based on the volume dependence of 
the statistical properties of the landscape |3q . 

Finally, within the Gaussian landscape model, it is pos- 
sible to exactly evaluate the Kauzmann curve Tr-(T, V) , 
the limit for the existence of the liquid. This curve is the 
locus of points where S con f vanishes, i.e., from Eq. J7J), 



{ 1 2a 2 (V) 
The following expression for Tk{T, V) results: 

T K (T,V) b 



(E -A)± y/2a 2 aN 



(8) 



(9) 



where the sign to be chosen is the one corresponding to 
the largest value solution of the equation. 



IV. MODEL AND SIMULATIONS 

We studied a system of 343 Lewis and Wahnstrom 
(LW) ortho-terphenyl model molecules, by means of 
molecular dynamics simulations in the (N, V, T) ensem- 
ble. The LW model is a rigid, three-site model with in- 
termolecular site-site interactions described by Lennard 
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FIG. 2: Volume dependence of the Gaussian landscape pa- 
rameters: a) a(V), b) E (V), and c) a 2 (V). 



Jones potential. The potential parameters are chosen to 
reproduce OTP properties such as its structure and dif- 
fusion coefficient [33. The integration time step for the 
simulation was 0.01 ps. With this model it is possible 
to reach very long simulation times; such long molecular 
dynamics trajectories allow us to equilibrate the system 
at temperatures below the temperature where the diffu- 
sion constant reaches values of order 10~ 10 cm 2 /s. We 
simulated 23 different densities for several temperatures, 
for an overall simulation time of order 10 fis. 

To calculate the inherent structures sampled by the 
system in equilibrium we perform conjugate gradient en- 
ergy minimizations to locate the closest local minima on 
the PEL, with a tolerance of 10 -15 kJ/mol. For each 
thermodynamical point we minimize at least 100 con- 
figurations, and we diagonalize the Hessian matrix of at 
least 50 configurations to determine the density of states. 
The Hessian is calculated choosing for each molecule the 
center of mass and the angles associated with rotations 
around the three principal inertia axis as coordinates. 

Further details on the numerical techniques used can 
be found in Ref . [13 ■ 
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FIG. 3: a) Inherent structure pressure, Pis- The IS pres- 
sure gets more negative on decreasing density, until a sudden 
jump upward takes place at around 225 cm 3 /mol. The jump 
signals that, during the minimization procedure (at constant 
volume), the maximum tensile strength has been overcome 
and a cavitation phenomenon has taken place, b) vibrational 
pressure, P V ib- Pvib is almost density independent. Only at 
densities close to the Sastry density, a V dependence is ob- 
served. 



POTENTIAL ENERGY LANDSCAPE 
PROPERTIES 



An analysis of the statistical properties of the land- 
scape for the LW ortho-terphenyl model has been re- 
cently performed in Ref . |17| . Here we expand such anal- 
ysis to lower and higher densities, with the aim of explor- 
ing the region of phase diagram where the steep repul- 
sive part of the potential is more relevant, and study the 
location of the liquid spinodal line and of the Sastry den- 
sity [2^, 0| . The larger density range considered allows 
us to estimate the volume dependence of the landscape 
parameters with great precision. 

The four panels of Fig. ^ show the results of the 
landscape analysis for two of the additional densities 
(198 and 200 cm 3 /mol), with the aim of confirming 
the possibility of describing the numerical data with 
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FIG. 4: Volume dependence of the total pressure P, inherent 
structures pressure Pis, and vibrational pressure P v ib for the 
isotherm T = 320 K. The total pressure P is monotonically 
decreasing, as expected for a system in equilibrium in the sta- 
ble liquid phase. The IS pressure shows instead a minimum 
around V = 222 cm 3 /mol, suggesting that, during the mini- 
mization process, a cavitation process has taken place in the 
V = 225 cm 3 /mol sample. As a result, the significant drop 
observed in P vi t is an artifact induced by cavitation. 



the Gaussian landscape model discussed above. Simi- 
lar data for other five densities can be found in Ref. 0] 
and are not shown here. Fig. ^a) shows the temper- 
ature dependence of the average inherent structure en- 
ergy (symbols) which, in agreement with Eq. @, can 
well be fitted by a 1/T law (solid line). Fig. □» 
shows the anharmonic entropy, evaluated according to 
a fit of the anharmonic contribution to the energy with 
a polynomial of third degree in T [13 . Fig. mc) 

shows the quantity ln(exp l^^ -3 ln(/3?iLjj(_E/s))|)' as 

a function of the basin depth energy Ejs- As pre- 
viously observed, an almost linear relation between 

ln(exp j^-f^M/^^/s))})' and E IS is found. To 
account for a minor curvature, a fit with a second order 
polynomial (Eq. [3J is reported. Finally, the Eis depen- 
dence of S con f is shown in Fig. dd). The parameters 
of the reported fit are constrained by the T— dependence 
of the parameters estimated for Eis |43l | . In agreement 
with Eq.s 10 and J7J), S con f is fitted by a second degree 
polynomial in Eis- 

As discussed in the previous section, from the T depen- 
dence of Eis and S con f , and from the Eis dependence of 

ln(exp |X)fcJi~ 3 m (/3^i(£'JS , ))|) / ; it is possible to evalu- 
ate the statistical properties of the landscape and their 
volume dependence, under the assumption of a Gaussian 
landscape. Fig. |2 shows the F-dependence of the param- 
eters a(V), E D (V), and & 2 (V). The parameter a shows a 
weak ^-dependence. From a theoretical point of view, we 
expect a to converge toward a constant value in the small 
volume limit, when the potential is essentially dominated 
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FIG. 5: Total system pressure P as a function of the volume 
V along the indicated isotherms. For each isotherm, the high- 
est volume shown is the volume where a jump in the inherent 
structure pressure is first observed, signaling bubble forma- 
tion -cavitation- in the system. The nucleation line (solid 
line) and the spinodal line (dashed line) are also shown. 

by the repulsive soft sphere r -12 part [l3j, and to increase 
on increasing V, to account for the larger configuration 
space volume. The F-dependence of E (V) is similar to 
the one found in other models; the presence of the min- 
imum is, indeed, connected to the progressive sampling, 
on compression, of the attractive part of the intramolec- 
ular potential, followed by the progressive probing of the 
repulsive part of the potential. As expected for simple 
liquids 50], the ^-dependence of cr 2 (V) is instead mono- 
tonic. 



VI. INHERENT STRUCTURES PRESSURE, Pis, 
AND VIBRATIONAL CONTRIBUTION, P vib 

Within the landscape approach, the pressure P can 
be exactly split in two contributions, one associated to 
the pressure Pjs experienced in the local minima (which 
are usually under tensile or compression stress), and a 
second contribution, P v ib, commonly named vibrational, 
even if a configurational part is also included as discussed 
in length in Ref. [U^l- Hence, in full generality, 

P{T, V) = P IS (T, V) + P vib (T, V). (10) 

Here, Pis(T,V) can be evaluated from the value of 
the virial expression in the inherent structure, while 
Pvib(T, V) can be evaluated as difference between P and 
Pis- 

Fig. Eta) shows the inherent structure pressure Pis 
and FigOJb) the vibrational pressure P V ih for several den- 
sities. The IS pressure gets more negative on decreas- 
ing density, until a sudden jump upward takes place at 
around 225 cm 3 /mol. The jump signals that, during the 
minimization procedure (at constant volume), the maxi- 
mum tensile strength has been overcome and a cavitation 
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FIG. 6: Inherent structures pressure, Pis, as a function of 
V, for several isotherms. Each curve has been shifted by 
n x 30 MPa, for a more clear presentation. Solid lines are 
the theoretical predictions according to the potential energy 
landscape equation of state of Ref. |3a| . The shadowed region 
marks the range of variability of Vsastry, the limit of stability 
for the inherent structures. Inset: Vsastry as a function of 
temperature for all the studied volumes is shown. Dashed 
line is the average value 224.9 cm 3 /mol. 



phenomenon has taken place. Therefore, the density at 
which the IS loses mechanical stability, recently named 
Sastry density (Vsastry)> is, f° r this model, close to 
224 cm 3 /mol. We also note that P vib (Fig. Efb)) is al- 
most density independent. Only at densities close to the 
Sastry density, a V dependence is observed. For densities 
lower than the Sastry density, both Pjg and P V ib do not 
reflect any longer bulk properties (being the local minima 
configuration affected by the presence of large voids). 

The previous observation stands out more clearly in 
Fig. 0] where the volume dependence of P, Pis and P„i6 
is shown for the isotherm T = 320 K. The total pressure 
P is monotonically decreasing, confirming that the stud- 
ied system is in the stable liquid phase up to V — 225 
cm 3 /mol. The IS pressure shows instead a minimum 
around V = 222 cm 3 /mol, suggesting that, during the 
minimization process, a small cavity has been created in 
the V — 225 cm 3 /mol sample. As a result, the significant 
drop observed at V = 225 cm 3 /mol in P v n, is an artifact 
induced by cavitation. It is also worth noting that a 
small decrease of P V ib is observed already at V = 222 
cm 3 /mol, suggesting that the cavitation phenomenon is 
preceded by a weak softening of the vibrational density 
of states on approaching the Sastry instability. 

VII. LIMITS OF STABILITY OF THE LIQUID 

We now focus on the limits of stability of the liq- 
uid state. Fig. [5] shows the volume dependence of the 
pressure for several of the studied isotherms. Cavitation 
marks the homogeneous nucleation limit for the system; 
it can be detected during the simulation by monitoring 




FIG. 7: Spinodal (open squares) and Kauzmann (open cir- 
cles) lines, together with their extrapolations (dashed lines) 
in the (T, V) plane. The relation between these two curves is 
discussed in the text. The vertical dashed line marks the Sas- 
try volume, the limit of stability for the inherent structures. 
Some iso-entropic lines are also shown. 



the time dependence of pressure and potential energy, 
which show a clear discontinuity when a gas bubble nu- 
cleates. After the cavitation, P increases and the po- 
tential energy decreases. We define the locus of homo- 
geneous nucleation as the largest volume, at fixed T, at 
which we managed to simulate the dynamics of a homo- 
geneous system. The fine grid of studied V values allows 
us to identify such locus with a significant precision (solid 
line). Although the calculated homogeneous nucleation 
line refers to a system composed of 343 molecules, it pro- 
vides an upper bound for larger systems. 

From the equilibrium P(T,V) data, i.e., in the range 
where no cavitation is observed, it is possible to estimate 
the volume at which P(V) has a minimum, by fitting the 
data according to the equation P m (V — V s ) 2 . In the 
mean field approximation, V s corresponds to the spinodal 
volume, and the temperature dependence of V s defines 
the spinodal locus (dashed line). 

We now focus on the information on liquid stability 
encoded in the inherent structures pressure Pis • Fig. 
shows, for several isotherms, Pis as a function of V. In 
analogy with the equilibrium data, a limit of stability for 
the inherent structures can be calculated estimating the 
minimum of Pjs ~ (V — Vsastry) are solid lines). 

It has been speculated that Vsastry represent the upper 
limit for glass formation, and that it may be identified 
with the T — > limit of the liquid-gas spinodal locus • 
Our calculations show that, as previously observed for 
other systems the Sastry volume does not depend 
significantly on the temperature (see inset in Fig. . 

We note that the Sastry volume is significantly smaller 
than the spinodal volume. The fact that Vsastry < V s 
implies a stabilizing role of the vibrational component 
of the pressure. Indeed, already in Fig. 0] it was shown 
that the stability region for P is larger than the one for 
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FIG. 8: Spinodal (open squares) and Kauzmann (open cir- 
cles) lines, together with their extrapolations (dashed lines) 
in the (P, T) plane. The vertical dashed line marks the Sastry 
volume. Some iso-entropic lines are also shown. 



Pis- This fact suggests that, close to the spinodal line, 
the vibrational component becomes volume dependent to 
compensate for the loss of stability arising from the Pis 
contribution. Clearly, if P vi b were V independent, then 
V s and Vsastry should coincide. There is also an interest- 
ing observation to make concerning the application of the 
potential energy landscape approach to liquids at large 
volumes. Indeed, between Vsastry and V s , the constant 
volume minimization procedure produces an inhomoge- 
neous IS configuration [E^]. When the inherent struc- 
ture contains voids, the determination of the landscape 
parameters proposed in this work becomes meaningless, 
and the link between the inherent structure and the cor- 
responding liquid state requires a more detailed model- 
ing. 

Next, we focus on the location of the relevant stabil- 
ity loci of the liquid phase under supercooling. Within 
the Gaussian landscape model, the limit of stability of 
the supercooled liquid is defined by the line at which 
the configurational entropy vanishes, the so-called Kauz- 
mann locus (Eq. jnj). As the spinodal curve is preempted 
by the homogeneous nucleation and, therefore, can not 
be approached in equilibrium, the Kauzmann locus can 
not be accessed due to the extremely slow structural re- 
laxation times close to it. Still, these two loci provide a 
characterization of the domain of stability of the liquid 
state, and retain some meaning in a limiting mean-field 
sense. 

Figs.[7|and[S]show the S con f = locus and the spinodal 
locus in the (T, V) and (P, T) planes. The (P, T) data 
have been calculated using the potential energy land- 
scape equation of state introduced in Ref. js^. Fig. 
also shows curves at constant configurational entropy, in 
the region where no extrapolations are required. In anal- 
ogy with the findings in Lennard Jones systems [29| , the 
volume where the two loci appear to meet is close to 
Vsastry, but at a temperature Tj ~ 73K different from 
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T = 0. Hence, the spinodal line terminates at a finite 
temperature, by merging with the S con f — line; this 
is at odds with what was suggested in Ref. |4flj |. For 
Tj ~ 73K, the glass will meet a mechanical instability on 
stretching; it is an interesting topic of research |29t l5.it l51| 
to understand the relations between the volume at which 
such instability takes place and Vsastry 

A recent thermodynamic analysis addresses the issue 
of the relative location of the two loci and the way these 
two loci intersect [3(1 . It has been suggested that the 
two lines must meet in (P, T) plane with the same slope. 
Making use of the equation of state of Ref. j3^ , the (T, V) 
data of Fig.[7|can be represented in the (P, T) plane. The 
phase diagram in (P, T) is shown in Fig. 03 where also 
Tk(P), the spinodal line and some iso-entropy curves are 
shown. An extrapolation of the low pressure behavior of 
the spinodal line shows that data are consistent with the 
possibility that the two lines meet with the same slope. 
If the pressure P(T) along the spinodal increases with T, 
as usually found in liquids, the meeting point defines the 
lowest temperature and pressure that can be reached by 
the liquid in equilibrium. These values are Tj ~ 73 K 
and Pr ~= -360 MPa. 

VIII. SUMMARY AND CONCLUSIONS 

In summary, we have studied the stability domain of 
the liquid state for a simple molecular model for ortho- 
tcrphcnyl. In particular, we have focused on the two 
limits of stability, one provided by the divergence of the 
structural relaxation times, the other one provided by 
the cavitation of the gas phase. Both of them have the- 
oretical mean-field bounds, the locus at which the con- 
figurational entropy vanishes and the locus at which the 
compressibility diverges, respectively. 

To evaluate the configurational entropy, we have de- 
veloped, along the lines of previous work for the same 
model fl7j , a potential energy landscape description of 
the system free energy, in the framework of the inher- 
ent structure thermodynamic formalism |18|. The land- 
scape analysis has required the evaluation of the statisti- 
cal properties of the landscape, which we have quantified 
in the volume dependence of total number, energy distri- 
bution, and relation between energy depth and shape of 
the potential energy landscape basins. 

The landscape analysis performed here is limited to 
volumes such that the minimization procedure, used for 
the evaluation of the inherent structures, results in an 



homogeneous structure. We have found that for volumes 
larger than the so-called Sastry density |5l| (Vsastry — 
225 cm 3 /mol) the system always cavitates upon min- 
imization. Cavitation prevents the possibility of esti- 
mating the statistical properties of the landscape for 
V > Vsastry This poses a serious problem to the ap- 
plication of landscape approaches in the version where 
minimizations are performed at constant V if the region 
close to the spinodal curve has to be investigated. In- 
deed, as shown in Fig. [7| the region of stability of the 
liquid phase extends well beyond Vsastry This suggests 
also that the region close to the spinodal curve is stabi- 
lized by vibrational contributions, which must overcome 
the destabilizing contribution arising from the configura- 
tional degrees of freedom (as discussed in Sec. IVIII . The 
stabilization in the vibrational properties appears to be 
accomplished by a softening of the vibrational density of 
states on approaching Vsastry We note on passing that, 
while in the (T, V) ensemble, estimates of the liquid free 
energy in the IS formalism are limited to V < Vs as try, 
formulations of the IS formalism in the (P, T) ensemble 
does not suffer from such limitation, since the minimiza- 
tion path would not meet any instability curve. 

In the framework of the IS formalism, we have esti- 
mated the locus at which configurational entropy van- 
ishes. The possibility of a finite T at which S con f van- 
ishes is encoded in the model selected to represent the 
data. For all densities and temperatures studied in the 
present work, the Gaussian landscape [2lJ well represents 
the data and provides a well defined Kauzmann locus, 
whose location in the phase diagram has been compared 
with the location of the spinodal line. We have shown 
that the spinodal and the S con f = loci may be extrap- 
olated to meet at Vsastry at a finite T. Data are also 
consistent with the possibility that, in the (P, T) plane, 
the two curves are tangent at Vsastry These two ob- 
servations are in agreement with the behavior recently 
predicted for hard and soft spheres complemented by a 
mean field attraction. It will be interesting to address 
in the future the relation between the field of stability 
of the liquid as compared to the field of stability of the 
glass state. 
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